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Abstract. We discuss results from numerical simulations of star cluster formation in the tur- 
bulent interstellar medium (ISM). The thermodynamic behavior of the star- forming gas plays 
a crucial role in fragmentation and determines the stellar mass function as well as the dynamic 
properties of the nascent stellar cluster. This holds for star formation in molecular clouds in 
the solar neighborhood as well as for the formation of the very first stars in the early universe. 
The thermodynamic state of the ISM is a result of the balance between heating and cooling 
processes, which in turn are determined by atomic and molecular physics and by chemical abun- 
dances. Features in the effective equation of state of the gas, such as a transition from a cooling 
to a heating regime, define a characteristic mass scale for fragmentation and so set the peak of 
the initial mass function of stars (IMF). As it is based on fundamental physical quantities and 
constants, this is an attractive approach to explain the apparent universality of the IMF in the 
solar neighborhood as well as the transition from purely primordial high-mass star formation to 
the more normal low-mass mode observed today. 
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of state - methods: numerical 



1. Introduction 

Identifying the physical processes that determine the masses of stars and their sta- 
tistical distribution, the initial mass function (IMF), is a fundamental problem in star- 
formation research. It is central to much of modern astrophysics, with implications rang- 
ing from cosmic re-ionisation and the formation of the first galaxies, over the evolution 
and structure of our own Milky Way, down to the build-up of planets and planetary 
systems. 

Near the Sun the number density of stars as a function of mass has a peak at a 
characteristic stellar mass of a few tenths of a solar mass, below which it declines 
steeply, and for masses above one solar mass it follows a power-law with an exponent 
dN/d\ogm oc m -13 . Within a radius of several kpc this distribution shows surprisingly 
little variation (Salpeter 1955; Scalo 1998; Kroupa 2001; Kroupa 2002; Chabrier 2003). 
This has prompted the suggestion that the distribution of stellar masses at birth is a 
truly universal function, which often is referred to as the Salpeter IMF, although note 
that the original Salpeter (1955) estimate was a pure power-law fit without characteristic 
mass scale. 

The initial conditions in star forming regions can vary considerably, even in the solar 
vicinity. If the IMF were to depend on the initial conditions, there would be no rea- 
son for it to be universal. Therefore a derivation of the characteristic stellar mass that 
is based on fundamental atomic and molecular physics would be highly desirable. In 
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this proceedings contribution we argue that indeed the thermodynamic properties of the 
star-forming cloud material determine the characteristic mass scale for fragmentation 
and subsequent stellar birth. The thermodynamic state of interstellar gas is a result of 
the balance between heating and cooling processes, which in turn are determined by fun- 
damental atomic and molecular physics and by chemical abundances. The derivation of 
a characteristic stellar mass can thus be based on quantities and constants that depend 
solely on the chemical abundances in a molecular cloud. It also explains why deviations 
from the "standard" mode of star formation are likely to occur under extreme envi- 
ronmental conditions such as those occurring in the early universe or in circum-nuclear 
starburst regions. 

2. Gravoturbulent Star Cluster Formation 

Stars and star clusters form through the interplay between self-gravity on the one hand 
and turbulence, magnetic fields, and thermal pressure on the other (for recent reviews 
see Larson 2003; Mac Low & Klessen 2004; Ballesteros-Paredes et al. 2006). Supersonic 
turbulence, even if it is strong enough to counterbalance gravity on global scales, will 
usually provoke local collapse. Turbulence establishes a complex network of interacting 
shocks, where converging shock fronts generate clumps of high density. These density 
enhancements can be large enough for the fluctuations to become gravitationally unstable 
and collapse, which can occur when the local Jeans length becomes smaller than the size of 
the fluctuation. However, the fluctuations in turbulent velocity fields are highly transient. 
The random flow that creates local density enhancements can disperse them again. For 
local collapse to actually result in the formation of stars, Jeans-unstable shock-generated 
density fluctuations must collapse to sufficiently high densities on time scales shorter 
than the typical time interval between two successive shock passages. Only then are they 
able to 'decouple' from the ambient flow and survive subsequent shock interactions. The 
shorter the time between shock passages, the less likely these fluctuations are to survive. 
Hence, the timescale and efficiency of protostellar core formation depend strongly on the 
wavelength and strength of the driving source as well as on the dynamic response of the 
gas, as defined by the equation of state, i.e. the balance between heating and cooling 
processes. 

The velocity field of long- wavelength turbulence is found to be dominated by large-scale 
shocks which are very efficient in sweeping up molecular cloud material, thus creating 
massive coherent structures. When a coherent region reaches the critical density for 
gravitational collapse, its mass typically exceeds the local Jeans limit by far. Inside the 
shock compressed region, the velocity dispersion is much smaller than in the ambient 
turbulent flow and the situation is similar to localized turbulent decay. These are the 
conditions for the formation of star clusters. The efficiency of turbulent fragmentation 
is reduced if the driving wavelength decreases. When energy is inserted mainly on small 
spatial scales, the network of interacting shocks is very tightly knit, and protostellar 
cores form independently of each other at random locations throughout the cloud and at 
random times. Individual shock-generated clumps have lower mass and the time interval 
between two shock passages through the same point in space is small. Collapsing cores 
are easily destroyed again and the resulting mass spectrum shows deviations from the 
observed IMF. All this points toward interstellar gas clouds being driven on large scales. 

Altogether, stellar birth is intimately linked to the dynamic behavior of the parental 
gas cloud, which governs when and where star formation sets in. The chemical and ther- 
modynamic properties of interstellar clouds play a key role in this process. In particular, 
the value of the polytropic exponent 7, when adopting an EOS of the form P oc p 7 , 
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strongly influences the compressibility of density condensations as well as the tempera- 
ture of the gas. The EOS thus determines the amount of clump fragmentation, and so 
directly relates to the IMF (Vazquez-Semadeni et al. 1996) with values of 7 larger than 
unity leading to little fragmentation and high mass cores (Li, Klessen, & Mac Low 2003; 
Jappsen et al. 2005). The stiffness of the EOS in turn depends strongly on the ambi- 
ent metallicity, density and infrared background radiation field produced by warm dust 
grains. The EOS thus varies considerably in different galactic environments (see Spaans 
& Silk 2000, 2005 for a detailed account). 



3. Formation of Stellar Clusters in the Solar Neighborhood 

Early studies of the balance between heating and cooling processes in collapsing clouds 
predicted temp eratures of the order of 10 K to 20 K, tending to be lower at the higher 
densities (e.g., Havashi fc Nakanol 19651 : Havashi Il966l : iLarsonl Il969l Il973bh . In their 
dynamical collapse calculations, these and other authors approximated this somewhat 
varying temperature by a simple constant value, usually taken to be 10 K. Nearly all 
subsequent studies of cloud collapse and fragmentation have used a similar isothermal 
approximation. However, this approximation is actually only a somewhat crude one, 
valid only to a factor of 2, since the temperature is predicted to vary by this much above 
and below the usually assumed constant value of 10 K . Given the str ong sensitivity of 



the results of fragmentation simulations like those of Li et al. ( 20031 ) to the assumed 



equation of state of the gas, temperature variations of this magnitude may be important 
for quantitative predictions of s tellar masses a nd the IMF. 

As can be seen in Fig. 2 of ILarsonl (|l985h . observational and theoretical studies of 
the thermal properties of collapsing clouds both indicate that at densities below about 
10 -18 g cm -3 , roughly corresponding to a number density of n = 2.5 x 10 5 cm -3 , the tem- 
perature generally decreases with increasing density. In this low-density regime, clouds 
are externally heated by cosmic rays or photoelectric heating, and they are cooled mainly 
by the collisional excitation of low-lying levels of C + ions and O atoms; the strong de- 
pendence of the cooling rate on density then yields an equilibrium tempe rature that 
decreases with increasing density. The work of iKovama fc Inutsukal (|2QQQh . which as- 
sumes that photoelectric heating dominates, rather than cosmic ray heating as had been 
assumed in earlier work, predicts a very similar trend of decreasing temperature with 
increasin g density at low den s ities. The three-dimensional magnetohydrodynamic simu- 
lations of iGlover fc Mac Low! (|2007l ) also produce a similar result, although in this case 
the point-to-point scatter is larger. The resulting temperature-density relation can be 
approximated by a power law with an exponent of about —0.275, which c orresponds to 
a polytropic equati on of state w ith 7 = 0.725. The observational results of Mversl (|l978l ) 
shown in Fig. 2 of lLarsonl (|l985h suggest temperatures rising again toward the high end 
of this low-density regime, but those measurements refer mainly to relatively massive 
and warm cloud cores and not to t he small, dense, cold cores in which low-mass stars 
form. As reviewed by lEvaneJ Jl999), the temperatures of these cores are typically only 
about 8.5 K at a density of 10 -19 gem -3 , consistent with a continuation of the decreasing 
trend noted above and with the continuing validity of a polytropic approximation with 
7 w 0.725 up to a density of at least 10 -19 gcm -3 . 

At higher densities, atomic line cooling becomes less effective as the cooling rates 
start to reach their local thermodynamic equilibrium (LTE) limits and as the line opaci- 
ties grow larger. Consequently, at densities above 10 -18 gcm -3 the gas becomes ther- 
mally coupled to the dust grains, which then control the temperature by their far- 
infrared thermal emission. In this high-density regime, dominated thermally by the dust, 
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there are few direct temperature measurements because the molecules normally observed 
freeze out onto the dust grains, but most of the available theoretical pr edictions are in 
good agreement concerning the expected thermal behavior of the g as ( Larson! Il973bl ; 
Low fc Lvnden-Belll Il976l : iMasunaga fc Inutsukal l2000l : iLarsonl [2005) . The balance be- 



tween compressional heating and thermal cooling by dust results in a temperature that 
increases slowly with increasing density, and the resulting temperature-density relation 
can be approximated by a power law with an exponent of about 0.075, which corresponds 
to 7 = 1.075. Between the low-density and the high-density regimes, the temperature is 
predicted to reach a minimum of 5 K at a density of about 2 x 10 -18 g cm -3 , at which point 
the Jeans mass is about 0.3 Mq. The actual minimum temperature reached is somewhat 
uncertain because observations have not yet confirmed the predicted very low values, but 
such cold gas would be very difficult to observe; various efforts to model the observations 
have suggested central temperatures between 6K and 10 K for the dens est observed 



prestellar cores, whose peak densities ma y approach 10 17 g cm 3 (e.g. IZucconi et al 
200 ll : lEvans et al.ll200ll : iTafalla et all feoO^. A power-law approximation to the equation 
of state with 7 « 1.075 is expected to remain valid up to a density of about 10 -13 g cm -3 , 
above which increasing opacity to the thermal emission from the dust causes the temper- 
ature to begin rising much more rapidly, resulting in an "opacity limit" on fragmentatio n 
that is somewhat below 0.01 M (|Low fc Lynden-B"ellll976l : lMasunaga fc Inutsukall200ol ). 

Adopting the piecewise polytropic equation of state outlined above, Jappsen et al. 
(2005) showed that the changing 7 from a value below unity to one somewhat above unity 
at a critical density n c influences the number of protostellar objects. If the critical density 
increases then more protostellar objects form but the mean mass decreases. Consequently, 
the peak of the resulting mass spectrum moves to lower masses with increasing critical 
density. This spectrum not only shows a pronounced peak but also a powerlaw tail 
towards higher masses. Its behavior is thus similar to the observed IMF. 

A simple scaling argument based on the Jeans mass Mj at the critical density n c leads 
to M c h oc n~ 0,95 . If there is a close relation between the average Jeans mass and the 
characteristic mass of a fragment, a similar relation should hold for the expected peak of 
the mass spectrum. The simulations by Jappsen et al. (2005) qualitatively support this 
hypothesis, but find a weaker density dependency M c h cx n~ 5±01 



4. Formation of Stellar Clusters in the Early Universe 

The formation of the first and second generations of stars in the early universe has far- 
reach i ng consequences for cosmic reioniza tion and galaxy formation (|Loeb fc Barkanal 
200 ll : iBromm fc Loebl I2QQ4J : Idoverl [2005). The physical processes that govern stellar 



birth in a metal- free or metal-poor environment, however, are still very poorly under- 
stood. Numerical simulations of the thermal and dynamical evolution of gas in pri- 
mordial protogalactic halos indicate that the metal-free first stars, the so called Pop- 
ulation II I, are expected to be very massive, with masses anywhere in the range 20 
2000 Ma (I Abel. Brvan. fc Nor man! 1 20021 : IBromm. Coppi. fc Larsonll2002l : lYoshida et ~ 
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the local IMF. 

This means that at some stage of cosmic evolution there must have been a transition 
from primordial, high-mass star formation to the "normal" mode of star formation that 
dominates today. The discovery of extr emely metal-poor subgiant stars in the Galactic 
halo with masses below one solar mass (jChristlieb et al.ll2002l : iBeers fc Christlieb 20051 ) 
indicates that this transition occurs at abundances considerably smaller than the solar 
value. At the extreme end, these stars have iron abundances less than 10 -5 Z , and car- 
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Figure 1. Three equations of state (E OSs) from Omukai et al. (2005) for different metallicities 
used bv lClark. Glover, fc Klessenl (|2007T ). The primordial case (solid line), Z = 10~ 6 Z© (dotted 
line), and Z = 10 -5 Z (dashed line), are shown alongside an example of a polytropic EOS with 
an effective 7 = 1.06. 



bon or oxygen abundances that are still < 10 ~ 3 the solar value. These stars are thus 
strongly iron deficient, which could be due to unusual abundance patt erns produced 
by enrichment from pair-instability supernovae ([Heger fc Wooslevl 120021) from Popula- 
tion III stars or due to mass transfer from a close binary companion ( Ryan et "all [2005; 
iKomiva et al. 2007 ). Ther e are hints for an inc reasing binary fraction with decreasing 
metallicity for these stars ( Lucatello et al]|2005f ). 

If metal enrichment is the key to the formation of low-mass stars, then logically there 
must be some critical metallicity Z cr i t at which the formation of low mass stars first 
becomes possible. However, the value of Z cr i t is a matter of ongoing debate. Some 
models suggest that low mass star formation becomes possible on ly once atomic fine- 
structure line cooling from carbon and oxygen becomes effective (iBromm et al] 12001 



Bromm fc Loebl 120031 : ISantoro fc Shulll 120061 : iFrebel. Johnson, fc Brommll2007l ). setting 
a value for Z cr i t at around 10 -3 5 Z . Another possibility is that low mass star forma- 
tion is a result of dust-induced fragmentatio n occurring at high densities, and thus at 
a very late stage in the protostellar collapse (ISchneider et al. 2002 : Omukai et al.ll2005l : 
Schneider et~al1l2006l : iTsuribe fc Omukaill2006h . In this model. 10~ 6 < Z crit < 10 -4 Z , 



where much of the uncertainty in the predicted value r esults from uncertainties in the 
d ust composition and the d e gree o f gas-phase depletion ([Schneider et al.1 [20021 . [2006). 



Clark. Glover, fc Klessenl ( 2007 ) modeled star formation in the central regions of low- 



mass halos at high redshift adopting an EOS similar to Omukai et al. (2005). They focused 
on a high-density regime with 10 5 cm -3 ^ n ^ 10 17 cm -3 . They find that enrichment 
of the gas to a metallicity of only Z = 10 -5 Z dramatically enhances fragmentation. A 
typical time evolution is illustrated in Figure O It shows several stages in the collapse 
process, spanning a time interval from shortly before the formation of the first protostar 
(as identified by the formation of a sink particle in the simulation) to 420 years afterwards. 
During the initial contraction, the cloud builds up a central core with a density of about 
n = 10 10 cm -3 . This core is supported by a combination of thermal pressure and rotation. 
Eventually, the core reaches high enough densities to go into free-fall collapse, and forms 
a single protostar. As more high angular momentum material falls to the center, the 
core evolves into a disk-like structure with density inhomogeneities caused by low levels 
of turbulence. As it grows in mass, its density increases. When dust-induced cooling 
sets in, it fragments heavily into a tightly packed protostellar cluster within only a few 
hundred years. One can see this behavior in particle density-position plots in Figure 
03 The simulation is stopped 420 years after the formation of the first stellar object 
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Figure 2. Time evolution of the density distribution in the innermost 400 AU of the proto- 
galactic halo shortly before and shortly after the formation of the first protostar at tsF- Only gas 
at densities above 10 10 cm -3 are plotted. The dynamical timescale at a density n — 10 13 cm -3 
is of the order of only 10 years. Dark dots indicate the location of protostars as identified by 
sink particles forming at n ^ 10 17 cm -3 . Note that without usage of sink particles to identify 
collapsed protostellar cores one would not have been able to follow the build-up of the proto- 
stellar cluster beyond the formation of the first object. There are 177 protostars when we stop 
the calculation at t = tsF + 420 yr. They occupy a region roughly a hundredth of the size of 
the initial cloud . With 18.7 Mq accreted at thi s stage, the stellar density is 2.25 x 10 9 M© pc -3 . 
Data are from lClark. Glover, fe Klessen (20071). 
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Figure 3. To illustrate the onset of the fragmentation process in the Z = 10 -5 Z© simulation, 
the graphs show the densities of the particles, plotted as a function of their x-position. Note that 
for each plot, the particle data has been centered on the region of interest. Results are plotted 
for three different output times, ranging from the time that the first star forms (£ s f) to 221 
years afterwards. The densities lying between the two horizontal das hed lines denote the range 
over w hich dust cooling lowers the gas temperature. The figure is from lClark. Glover. &; Kle ssen 
(f2007h . 
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Figure 4. Mass functions resulting from simulations with metallicities Z = 10 -5 Z© (left-hand 
panel), Z = 10 -6 Z© (center panel), and Z = (right-hand panel). The plots refer to the point in 
each simulation at which 19 M© of material has been accreted (which occurs at a slightly different 
time in each simulation). The mass resolutions are 0.002 M© and 0.025 M© for the high and low 
resolution simulations, respectively. Note the similarity between the results of the low-resolution 
and high-resolution simulations. The onset of dust-cooling in the Z = 10 -5 Z© cloud results in 
a stellar cluster which has a mass function similar to that for present day stars, in that the 
majority of the mass resides in the lower-mass objects. This contrasts with the Z = 10 -6 Z© and 
primordial clouds, in which the bulk of the cluster mass is in high- mass stars. The figure is from 
iClark. Glover, fe Klessen] (p007h . 



(sink particle). At this point, the core has formed 177 stars. The evolution in the low- 
resolution simulation is very similar. The time between the formation of the first and 
second protostars is roughly 23 years, which is two orders of magnitude higher than the 
free-fall time at the density where the sinks are formed. Note that without the inclusion 
of sink particles, one would only have been able to capture the formation of the first 
collapsing object which forms the first protostar: the formation of the accompanying 
cluster would have been missed entirely. 

The fragmentation of low-metallicity gas in this model is the result of two key features 
in its thermal evolution. First, the rise in the EOS curve between densities 10 9 cm -3 
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and 10 11 cm -3 causes material to loiter at this point in the gravitational contraction. A 
similar behavior at densities around n = 10 3 cm -3 is discussed by Bromm et al. (2001), 
who call it a loitering phase. The rotationally stabilized disk-like structure, as seen in 
the plateau at n « 10 10 cm -3 in Figure [3l is able to accumulate a significant amount 
of mass in this phase and only slowly increases in density. Second, once the density 
exceeds n ~ 10 12 cm -3 , the sudden drop in the EOS curve lowers the critical mass for 
gravitational collapse by two orders of magnitude. The Jeans mass in the gas at this stage 
is only Mj = 0.01 M . The disk-like structure suddenly becomes highly unstable against 
gravitational collapse and fragments vigorously on timescales of several hundred years. A 
very dense cluster of embedded low-mass protostars builds up, and the protostars grow 
in mass by accretion from the available gas reservoir. The number of protostars formed 
by the end of the simulation is nearly two orders of magnitude larger than the initial 
number of Jeans masses in the cloud set-up. 

Because the evolutionary timescale of the system is extremely short - the free-fall time 
at a density of n = 10 13 cm -3 is of the order of 10 years - none of the protostars that 
have formed by the time that the simulation is stopped have yet commenced hydrogen 
burning. This justifies neglecting the effects of protostellar feedback in this study. Heating 
of the dust due to the sign ificant accretion luminosities of the newly-formed protostars 
will occur (|Krumholdl2006h . but is unlikely to be important, as the temperature of the 
dust at the onset of dust-induced cooling is much higher than in a typical Galactic 
protostellar core (Td us t ~ 100 K or more, compared to ~ 10 K in the Galactic case). The 
rapid collapse and fragmentation of the ga s also leaves no time for dynamo amplification 
of magnetic fields (|Tan fc Blackmanll2004f h which in any case are expec ted to be wea k 
and dynamically unimportant in primordial and very low metallicity gas (|Widrowll2002l ). 

The forming cluster represents a very extreme analo gue of the clustered s tar formation 
that we know dominates in the present-day Universe rtLada fc Ladall2003h . A mere 420 
years after the formation of the first object, the cluster has formed 177 stars (see Figure 
[2|). These occupy a region of only around 400 AU, or 2 x 10 -3 pc, in size, roughly a hun- 
dredth of the size of the initial cloud. With ~ 19 M accreted at this stage, the stellar 
density is 2.25 x 10 9 M pc -3 . This is about fi ve orders of magnitude greater th an the stel- 
lar density in the Trapezium cluster in Orion ( Hillenbrand fc Hartmannlll998l ) and about 
a thou sand times greater than that in the core of 30 Doradus in the Large Magellanic 
Cloud ([Massev fc Hunterlll998h . This means that dynamical encounters will be extremely 
important during the formation of the first star cluster. The violent environment causes 
stars to be thrown out of the denser regions of the cluster, slowing down their accretion. 
The stellar mass spectrum thu s depends on both the details of the initia l fragmentation 
process (e.g. as discussed bv ljappsen et al. 2005 : Clark fc Bonnelll 2005|) as well as dy - 
namical effects in the growing cluster (Bo nnell et al. ! 20011: iBonnell. Bate fc Vinel koCte). 
This is different to present-day star formation, where the situation is less clear-cut and 
the relative importance of these two processes may vary strongly from region to region 
dKrumholz. McKee. fc Kleinl 120051 : IBonnell fc Bate! 120061 : IBonnell. Larson fc Zinnecker ' 
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The mass functions of the protostars at the end of the Z = 10 -5 Z simulations 
(both high and low resolution cases) are shown in Figure [H (left-hand panel). When 
the simulation is terminated, collapsed cores hold ~ 19 M of gas in total. The mass 
function peaks somewhere below 0.1 M and ranges from below 0.01 M to about 5 M . 
This is not the final protostellar mass function. The continuing accretion of gas by the 
cluster will alter the mass function, as will mergers between the newly-formed protostars 
(which cannot be followed using our current sink particle implementation). Protostellar 
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feedback in the form of winds, jets and Hll regions may also play a role in determining 
the shape of the final stellar mass function. However, a key point to note is that the 
chaotic evolution of a bound system such as this cluster ensures that a wide spread of 
stellar masses will persist. Some stars will enjoy favourable accretion at the expense of 
others that will be thrown out of the system (as can be seen in Figure [2]), thus having 
their accretion effectively terminated (see for example, the discussions in Bonnell & Bate 
2006 and Bonnell, Larson & Zinnecker 2007). The survival of some of the low mass stars 
formed in the cluster is therefore inevitable. 

In the Z = and Z 10~ 6 Z calculations IClark. Glover, fc Klessenl (|2QQ7f ) find that 
fragmentation of the gas occurs as well, albeit at a much lower level than in the Z = 
1O _5 Z0 run. The mass functions from these simulations are shown in Figured! (middle and 
right-hand panels), and are again taken when ~ 19 M of gas has been accreted onto the 
sink particles, the same amount as is accreted by the end of the Z = 1O _5 Z0 calculations. 
Both distributions are considerably flatter than the present day IMF, in agreement with 
the suggestion that Population III stars are typically very massive. The fragmentation 
in the Z = 10 -6 Z simulation is slightly more efficient than in the primordial case, with 
33 objects forming. 



5. Conclusions 

In this proceedings paper we have discussed several studies, where the thermodynamic 
behavior of the interstellar medium plays a crucial role in fragmentation and subsequent 
stellar birth. These examples range from star formation in the solar neighborhood at the 
present day to the formation of the very first stars in the early universe, and support 
the idea that the distribution of stellar masses depends, at least in part, on the thermo- 
dynamic state of the star-forming gas. Dips in the effective EOS, such that the relation 
between temperature T and density p changes from decreasing T with increasing p to 
increasing T with increasing p, i.e. the transition from a cooling (7 < 1) to a heating 
(7 > 1) regime, define a characteristic mass scale for fragmentation. 

The thermodynamic state of interstellar gas is a result of the balance between heating 
and cooling processes, which in turn are determined by fundamental atomic and molecular 
physics and by chemical abundances. The derivation of a characteristic stellar mass can 
thus be based on quantities and constants that depend solely on the chemical abundances 
of the star forming gas. This is an attractive feature explaining the apparent universality 
of the IMF in the solar neighborhood as well as the transition from purely primordial 
high-mass star formation to the more normal low-mass mode observed today. Clearly 
more work needs to be done to investigate the validity of this hypothesis. 
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